Late Time Analysis for Maximal Slicing of Reissner-Nordstrom Puncture Evolutions 
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We perform an analytic late time analysis for maximal slicing of the Reissner-Nordstrom black 
hole spacetime. In particular, we discuss the collapse of the lapse in terms of its late time behavior 
at the throat and at the event horizon for the even and the puncture lapse. In the latter case we 
also determine the value of the lapse at the puncture. Furthermore, in the limit of late times slice 
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numerical experiments for a Schwarzschild black hole with puncture lapse and find agreement with 
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I. INTRODUCTION 

Choosing a lapse function that results in maximal slic- 
ing is a popular coordinate choice when decomposing the 
4-dimensional Einstein equations into (3-|-l)-dimensional 
form. This coordinate choice corresponds to the condi- 
tion that the mean extrinsic curvature of the slices van- 
ishes at all times, leading to a maximal volume property 
of the hypersurfaces 

For a single Schwarzschild black hole, maximal slices 
can be constructed analytically 0, Q ■ One can obtain 
spacelike hypersurfaces which extend from the "right- 
hand" spatial infinity to the "left-hand" inner infinity of 
the extended Schwarzschild spacetime. These slices give 
a complete foliation outside the event horizon while ap- 
proaching a limiting slice inside that does not reach the 
physical singularity. In this sense maximal slices avoid 
the physical singularity, and even in more general situ- 
ations maximal slices have been found to be singularity 
avoiding 0|. 

For these reasons, maximal slicing has been used fre- 
uently in numerical simulations of one black hole, e.g. 
1^ [3, and also in binary black hole mergers in full 
3D, e.g. B, B . Furthermore, there now are shift con- 
ditions [lOj which in many cases overcome the so- 
called slice stretching problem that previously limited 
black hole evolutions like those in to short evolution 
times, although we do not consider such shift conditions 
here. 

In a previous paper |lll | we established the connection 
between the analytically known solutions for maximal 
slicing of the Schwarzschild spacetime 0,0 and the max- 
imal slicing computed in the puncture evolution method 
for black holes 0. Only "odd" and "even" boundary 
conditions for the lapse had been considered before, but 
the numerical lapse of the puncture evolution is not of 
this type and in the past it was not clear whether the 
numerically obtained maximal slices exist analytically. 
We proved by explicit construction that there exists a 
maximal slicing of the Schwarzschild spacetime such that 



the lapse has zero gradient at the puncture, which is the 
boundary condition that has been observed to hold in 
numerical evolutions. We referred to this lapse as punc- 
ture lapse or "zgp" lapse and presented our results for 
the non-extremal Reissner-Nordstrom metric, generaliz- 
ing previous constructions of maximal slices. 

However, the analytic expressions we found in |llj are 
not very suitable for a discussion of the late time be- 
havior of the maximal slices. In the present paper, in 
order to provide information about the 4-metric in the 
limit of late times, we perform a late time analysis (LTA) 
for maximal slicing with even orpuncture lapse. This 
study is based on the results in [ill lia | using a method 
pointed out by R. Beig and N. 6 Murchadha Q. In the 
latter reference the late time behavior of the even lapse at 
the throat was derived analytically for the Schwarzschild 
metric. Showing that the lapse collapses exponentially 
there on a time scale given by = |-\/6M « 1.8371M, 
their analytical result confirmed earlier numerical results 
by L. Smarr and J.W. Yorkjll. 

Our LTA will extend |3| by introducing electrical 
charge to cover the Reissner-Nordstrom spacetime. Fur- 
thermore, the analysis will be extended to four locations 
of interest in numerical relativity, namely "from left to 
right" - to be understood in terms of a Carter-Penrose 
diagram or the numerical grid (c.f. Fig. 4 or Fig. 5 in [TH 'l 
- the puncture, the left-hand event horizon, the throat, 
and the right-hand event horizon. Studying the behavior 
of even and puncture lapse at those "markers" is a useful 
first step in the discussion of the "collapse of the lapse" . 

In addition, we will study the 3-metric in the limit of 
late times. Here, by deriving the late time behavior of the 
event horizon in volume radius coordinates, we are able to 
show that a conjecture made in Q for even boundary con- 
ditions does hold there. Furthermore, for isotropic grid 
coordinates used in puncture evolutions with vanishing 
shift, we determine slice stretching effects [sj. Whereas 
those effects frequently have been observed numerically 
[^^3, to our knowledge they are investigated here on an 
analytic level for the first time. We discuss the develop- 
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mcnt of a rapidly growing peak in the radial component of 
the metric ("slice wrapping") and the "outward" -drifting 
of coordinate locations such as the right-hand event hori- 
zon ("slice sucking"). 

The paper is organized as follows. In Sec. ^ we recall 
the analytic solution for maximal slicing of the Reissner- 
Nordstrom spacetimc and discuss briefly odd, even, and 
"zgp" boundary conditions. In Sec. IIIII we derive the 
behavior of the even and puncture lapse at late times an- 
alytically and compare with numerical results. In Sec. lIVI 
slice sucking taking place at the event horizon is studied 
for volume radius coordinates. In addition, for isotropic 
grid coordinates we discuss the behavior of the 3-metric 
in the limit of late times in order to point out slice stretch- 
ing effects. Again we compare with numerical results. We 
conclude in Sec. The expansions underlying our LTA 
arc proven in the Appendix. 



radial component of the 3-metric can be expressed as 

/3(r,r) = ^^^Ii^^C(r) (5) 

and 

1 

/(r) + ^ Pc{T,r) 

respectively. 

Next we want to re-derive the odd and even lapse and 
discuss a numerically motivated lapse with zero gradient 
at the puncture, the "zgp" or puncture lapse. 



B. Odd Lapse 



II. MAXIMAL SLICING OF THE 
REISSNER-NORDSTROM SPACETIME 

A. General solution in radial gauge 

In this Section we briefly summarize the formulas 
stated in [Tll | for maximal slicing of the Reissner- 
Nordstrom spacetime. Following Estabrook et al. we 
use the radial gauge, 

ds^ = {-a^ + —)dT'^ + 2l3dTdr + -idr'^ dn^, (1) 
7 

to solve the Einstein-Maxwell equations in spherical sym- 
metry under the condition of maximal slicing being van- 
ishing trace of the extrinsic curvature, K = jijK^^ = 0. 
We showed that the general solution of the lapse satis- 
fying the corresponding elliptic maximality condition is 
given by 

a[T,r) = 2 + -r / — ■ (2) 

V dr J,. pc{T,y)2j 

Here C and D are functions of r only and the polynomial 

pcir^r) = r^f{r) + C^t) = - 2Mr^ + qM\^ + C^{t) 

(3) 

based on 

has been introduced. As in we consider the non- 
extremal Reissner-Nordstrom spacetime. So M > de- 
notes the mass of the black hole and the dimension- 
less charge parameter defined hy q = is restricted to 
< q < 1. 

Imposing boundary conditions on the lapse fixes the 
so far undetermined functions C and D. However, inde- 
pendently of these boundary conditions the shift and the 



The first particular lapse choice is referred to as odd 
(in some references also as antisymmetric) because of its 
underlying antisymmetry with respect to the throat. It 
corresponds to the boundary conditions 

1™ a odd = ±1 ^'^odd- (7) 

r — ^oo 

To measure proper time at infinity (on the right-hand 
side of the throat denoted by a superscript "-f " ) the lapse 
is unity there and due to the antisymmetry the lapse 
is minus one at the puncture (i.e. at the compactified 
infinity on the left-hand side denoted by "— "). The odd 
lapse, 

■'i.AC,r)=±f^^±^. (8) 

is obtained by setting in Q D{Todd) ~ ±1 and ~ 0. 
Note that the slice label C can be chosen indepenciently 
of the time at infinity Todd- Furthermore, whereas event 
and Cauchy horizon with r± = (1 ± ^/l — q)M are given 
by the real roots of /(r) and are independent of C, the 
throat rc is obtained as a real root oipci''') and therefore 
depends on the choice of C . In particular, this implies 
C = when starting with initial data where the throat 
coincides with the event horizon. In this case we obtain 
from © the so-called Schwarzschild lapse. 



"^dd(^^ = o-0 = ±V7M, (9) 

from |SJ) a vanishing shift, and from (O the radial met- 
ric component 7(r) ~ Hence the static Reissner- 
Nordstrom spacetime written in Schwarzschild coordi- 
nates, 

ds^ = -f{r) df + -i- dr^ + r^dn^, (10) 
fir) 

is recovered. 
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Applying for r > r_|- a purely spatial coordinate trans- 
formation, 



R 



R + M 



M 
2R 



4i?2 

(1 - qW^ 
iR 



(11) 



chosen such that r and R coincide at infinity and r_(_ is 
mapped to i?+ = ^y/T~~qM , the isotropic lapse 

a„,,(C - 0, i?) - (2i? + M)2_^Ap (12) 

is found. The shift, of course, still vanishes and the 
3-metric is given by '^^Us^ = ^^{R){dR^ + R^ dVi^) with 
conformal factor 



M 
2R 



(13) 



One can readily verify that for isotropic coordinates the 
isometry 



R 



(1 ~ q)M 
AR 



/f2 



(14) 



is present. In particular, i? = is simply a compactified 
image of infinity which is referred to as puncture. 



where the integral 

Kc{r) = 



y{y-MI+^-^)dy 



(17) 



is defined for r > rc- Here the even height function is 
given by 



teveniC,r) = Hc{r) 



with the integral 



r 

Hc{r) = - j 



Cdy 



f{y)Vpc{y) 



(18) 



(19) 



again defined for r > rc- Note that in H19|l the integra- 
tion across the pole at r_|_ is taken in the sense of the prin- 
cipal value and the corresponding slices extend smoothly 
through the event horizon r+ and the throat rc- Since 
proper time is measured at spatial infinity, from 118|l in 
the limit r ^ oo one can infer 



i(C) lim t^yeniC^r) ^ Hdoo)- 



(20) 



In addition, we want to mention that in deriving these 
expressions (see 0| for details) the statement 



C. Even Lapse 

Our next example is the even (in some references also 
called symmetric) lapse satisfying the boundary condi- 
tions 



lim Q:J^,e„ = 1 yTeven- 



(15) 



Here in (jSJ D{Teven) is fixed to be unity and as pointed 
out in ^] the time-dependence manifest in C(Te„e„) can 
be determined by imposing the requirement of smooth- 
ness across the Einstein-Rosen bridge when passing 
to the line element (IIUII of the Reissner- Nordstrom space- 
time in Schwarzschild coordinates. 

In [Tll | we have derived expressions for the even lapse 
and the corresponding "height function" which generalize 
the formulas stated for the Schwarzschild case in Q by 
the electric charge. It turns out that the even lapse can 
be written as 



dC dU 



dTp, 



dC 
1 



(16) 



^Kcir)] , 



A"-^-^^ 2(.-iM;^)v^ 4^-^-(-) 

has been used, which in the limit r ^ oo yields the equa- 
tion 



^Hc(oo) = -^Kcioo) 



(22) 



to be used later on. Furthermore, we want to point out 
that as Teven IS running from zero to infinity, the slice 
label C takes values ranging from zero to Cum ■ As shown 
in [1^] the latter is determined by pc{tc) = for 

rcu^ = 4(3 + v/9^)Af '-^ — (23) 
to be given by 



^ a/27 - 36g + 8g2 + (9 - d>qf/-^M'^ 



^ -V3Af2 _ i.2990M^ 



(24) 



With rcii^ > the throat rc - being the "innermost" 
point on a slice labeled by C - never reaches the singu- 
larity at 7' = 0. Hence the singularity avoidance of maxi- 
mal slices for the Rcissner-Nordstrom spacetimc becomes 
apparent, c.f. corollary 3.7 of |^. 
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D. Lapse with zero gradient at the puncture 

Next we want to discuss a numerically motivated 
boundary condition arising for puncture data. Here a 
"radial grid coordinate" x measures the Euclidean dis- 
tance from the puncture situated at the origin of a spa- 
tial grid and r = x for {r — > oo, x — > 00} and r oc ^ for 
{r — > 00, a; 0} hold. On the initial slice this grid co- 
ordinate is identical to the isotropic radial coordinate R 
as in Hll() . but x differs at later times depending on the 
boundary conditions and the shift used in the numerical 
simulation. Demanding as before unit lapse at right-hand 
infinity, 



lim a"*" 



lim 



1 Vr, 



(25) 



one can label the slices with proper time. Numerically, 
the inner boundary condition arises as a special solution 
to the maximal slicing equation in the presence of a coor- 
dinate singularity in the metric. No boundary is actively 
enforced, but a numerical algorithm finds a solution with 
sufficiently fast fall-off in the derivative of the lapse such 
that, even though there is a coordinate singularity in the 
metric, the lapse is regular there. In particular, for the 
lapse a vanishing gradient at the puncture is obtained 
which can be formulated as 



As pointed out in the time-dependence is consis- 
tently derived from 



d 



dC ^^'P 



(30) 



by the "zgp" height function 



t%^{C, r) = Heir) ± {Hc{^) - ^) (31) 



when measuring time at right-hand spatial infinity by 
TzgpiC) = ^Inn i+p(C,r) = 2i/c(oo) - ^. (32) 

But note that the time measured at the puncture ap- 
proaches the finite value 



lim lim t JC,r) 



lim- 



M 



(33) 



III. 



LATE TIME BEHAVIOR OF THE EVEN 
AND PUNCTURE LAPSE 



lim dxOi^n-n oc lim r^dra^„„ ~ Vr^ 

x-fO r->oo 



(26) 



making use of dx oc r^dr in this limit. 

Since the maximality condition determining the lapse 
is linear in a, one can derive as in [Tll | the puncture lapse 
as a superposition, 

a± p(C, r) = $(C) . a± ,„(C, r) + (1 - ^C)) ■ a± ,(C, r), 

(27) 

of the linearly independent odd and even lapse. Impos- 
ing the boundary condition ]2<il , the "multiplicator func- 
tion" <!> present in (|27|l is given by 



1 Kcjoo) 

2 Ac (00) 



(28) 



M 



As we will show later, the integral Kc{oo) diverges in the 
limit of late times and $ hence approaches the value i. 
For this reason one can think of the "zgp" lapse as being 
essentially an average of the odd and even lapse. 

Combining the results ©, l(T^ . and we can ex- 
press the "zgp" lapse as 



± N VpcF) dC dt, 
otznp{C^r) = - 



(29) 



2Ac(oo)-f ^ V,_|M+S|li 



^(7.c(0±(^c(oo) + -|)) 



A. Motivation for a late time analysis 

In practice one is interested in statements for the 
4-metric (and in particular for the lapse profile) at a given 
time at infinity f. This time can be chosen without loss 
of generality to be < f < cx) since for f < one could 
redefine the sign of the time coordinate and for f = 
and f = 00 the solutions corresponding to C = and 
C = Ciim are known analytically. So as a first step one 
has to invert the relation t(C) - stated for even boundary 
conditions in (|20|1 and for "zgp" boundary conditions in 
()32|l - in order to find the appropriate slice label C with 
which our analytic expressions are to be evaluated in a 
second step. Both steps, however, arc non-trivial since 
the primitive for the integrals Hc{r), H19(l . and Kc{r), 
((T7|l . can not be written in closed form. Even worse, these 
integrals are numerically ill-conditioned in the sense that 
they diverge in the "late time limit" C Cum , as the 
Schwarzschild radius at the throat rc becomes with re, 
a double counting root of the polynomial pc (r) , the root 
of which appears in the denominator of the integrands 
of Hc{r) and Kci^)- An efficient numerical implemen- 
tation to calculate these integrals following Appendix 3 
of |T3| can be found in but the algorithm described 
there breaks down for C being too close to either or 

The problem of calculating for a given "late" time the 
appropriate slice label is apparent from Fig. ^ Concen- 
trating in particular on the "zgp" Schwarzschild case, we 
can show (using results of the now following LTA) that 
the sfice labels °C for r,™ = lOM, lOOM and lOOOM 



5 



differ from '^Cum 
by only 



r\/3M^ corresponding to r^gp = oo 



lim 



Czgp{Tzgp 
^Czgp{f: 



lOM) < 9 • 10-'"' 
lOOM) < 5 • 10 
lOOOM) < 9 • 10-^^'' 



24 M\ 



(34) 



Due to these difficulties when evaluating the analytic 
expressions, hypersurfaces corresponding to f > lOA/ 
are rarely shown in spacetime or embedding diagrams 
found in the literature for the maximallysliced Reissner- 
Nordstrom or Schwarzschild spacetime 0, H, 0, 0| ■ 



TIM] 



0.2 0.4 




C[Ml 



FIG. 1: On a logaritlimic scale the time at infinity corre- 
sponding to even and "zero gradient at the puncture" bound- 
ary conditions, r^vm (gray color) and r^gp (black), is shown as 
a function of the slice label C for the four values of the charge 
parameter q = ^ being (solid), 0.5 (dashed), 0.99 (dot- 
dashed), and 1 (dotted). Note that the range < C < ''Cum 
is getting smaller as ''Cum declines from its Schwarzschild 
value °Ciim = jVSM^ with increasing q such that ^Ciim — 0. 
The steep gradient in t(C) present as C ^ Cum causes nu- 
merical problems in finding C for a given time at infinity r as 
described in the text. 



B. Concept of a late time analysis 

For practical purposes one is hence interested in ap- 
proximations for the C-dependence of r. The main idea 
of any LTA now is to introduce an additional time coor- 
dinate which is particularly suitable when discussing the 
limit of late times. 

It is interesting to point out that in Eq. (A3. 16a- c) of 
[l7j for the even Schwarzschild case discussed there, an 
approximately linear relationship 



-0.918555 In 



lim 



c 



AP 



M + 0.2865Af (35) 



between the time at infinity Teven and the "logarithmic 
Estabrook-Wahlquist time" defined by - In M 
was found empirically. Anticipating again results 



Pi -0.918559M and 



of our LTA, this statement turns out to be a 
very precise approximation for "late times" where 
the constants appearing in H35() can be found an- 
alytically to be given by 

°A-H(is^ - In [M]fn « 0.286424M. 

However, we will not follow up on this logarithmic 
Estabrook-Wahlquist time further, but instead adopt the 
approach of Q , which can be summarized as follows. Our 
LTA is based on 5 introduced by 



5 



rc ~ rci. 



(36) 



as the difference between the radial coordinate at the 
throat rc on a slice labeled by C and its final value 
'"Clin, corresponding to C^m- Obviously, S is lying in 
< (5 < (7+ when defining 



<t7± = r± 



M 

3 J?/ 
2 



(37) 



and S approaches zero in the limit of late times. 

By expanding the proper time at infinity - which we 
have found for even and "zgp" boundary conditions in 
(|20|1 and H32|) - and inverting the relationship, an ex- 
ponential decay of S with r is found. For this decay a 
"fundamental time scale" turns out to be 



o<n = - 



Cii 



1 M- 



2 ^-aj^a^ ^ 

1.8371M, 



4 



(38) 



where is defined in terms of the following quantities: 
t = (1 + v/9~~8^)A/ ^ 



< 
< 



0^ 



%q + 3v/9 - 8g)A/ 
2-q 



2 9^0 9 



o< ?d 



A± 



2AI 




(39) 
(40) 

(41) 



In 



±- 



- o-±^ + 1 
a±^ + 2(i^ - 



+ CT±/1 + I') 



<7±{n + 2(cr± + Wcr± + (T±H 



q^O 



X± 



3*3 '» 







9V6+22 
5+3^3 



M 



2.2528ji- 



2q± 





9(-3 + 2g) 
V9~~8^ 



^)) 
(42) 

(43) 



Note that a time scale ^Vl w 1.82M has been found nu- 
merically for the Schwarzschild spacetime in [l^ and that 



6 



1.8371M has been de- 



its analytical value ~ |\/6Af 
rived in 

Then if one is interested in studying the lapse, one has 
to expand the lapse functions, Hlt)|) or H29|l . in terms of 
S and discuss the limit (5 — > 0. When calculating these 
expansions the following important expressions arise: 



> 77 



In 



{M 7^)2 



(T+(/i + 2^) 



16(2cr_ - i\/ + 2y'cr2 + cr+^ + zy) 
(2(T_ -il/)2(M-^)2 



3cr|/i + {2v - 3(7+ /.i) 4/0-2 _|_ ^^^^ _|_ 



1 7V21n 



36 



(2a+t.)2 

9%/6 + 22 
3V2 + 4 



-0.7385 



6(3^3 + 2) j ^ 
(44) 



and 



A 



-r21n 



2^ + H 
8^^ 



17 In [M] (-Vein 18(3^2-4) 



-2 In 



3V3 - 5 



9^6-22 
1.8371il/ln[M] -0.2181M. 



)M 



(45) 



Here A is playing the role of the constant A as defined in 
Eq. (3.41) of 0. As the reader might guess by now, the 
calculations are rather involved. For this reason we have 
only performed the (5-expansions at the puncture, the 
left-hand event horizon, the throat and the right-hand 
event horizon. But information obtained at these mark- 
ers essentially allows us to sketch the lapse profile and 
discuss in leading order the "collapse of the lapse". We 
have summarized our results for even and "zgp" bound- 
ary conditions in Subsecs. IIII CI and IIII Dl respectively. 
The proof of those statements can be found in the Ap- 
pendix. 



C. Late time behavior of tlie even lapse 

Expanding proper time at infinity in 5 and inverting 
the relation, 5 is found to decay exponentially with time 

as 



5^ 



, exp 



^eve 



©(exp 



-2- 



^eve 



)• (46) 



We would like to refer to the time scale 

Vleven = n -VgM « 1.8371M (47) 



defined in H38|) as being the fundamental time scale 
which characterizes the collapse of the lapse when maxi- 
mally slicing a Rcissner-Nordstrom black hole. The pre- 
exponential factor in l|46|l turns out to be 



0.8880M. (48) 



Due to its boundary conditions the even lapse is one at 
both puncture and infinity. Expanding the even lapse in 
d at the throat, which is the center of symmetry, yields 





'A' 




"OA" 


^even — exp 


n 


— > exp 


on 



2V2 S 



3 M 



0.9428- 



M 



(49) 



We want to mention here that for the Schwarzschild 
spacetime this decay of 0{6) at the throat has been de- 
rived already in Q . Extending the latter study we obtain 
for the lapse at the event horizon the expansion 



G 



r 



+ 

M - fi 



X- 



~2" (GlimT^ 



) 5^ + 



16 



V3- 



1 



:( 7\/6 



In 



-330 



0.3248+ 1.2265 



144^2 V 
12\/3- 36\/6 
52 



9V6 + 22 



3\/2 

il/2 



M2 



4 

0{8^), (50) 



which because of symmetry is identical for the left-hand 
and right-hand event horizon. Note that the finite value 
approached there as ymp totically is in agreement with our 
conjecture made in |ll| . 



D. Late time behavior of the puncture lapse 



Similarly, the late time behavior of 5 as a function of 



T2„„ is determined by 



^zgp exp 



gp 



©(exp 



-2 



zgp 



), (51) 



where the time scale of the exponential decay is found to 
be twice the fundamental time scale, 



2n ^ -VQM w 3.6742Af. (52) 
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Here the prefactor Kzgp can be calculated as 



exp 



exp 



A — In J- 

2 A'/ 



OA 



on 



0.6236M. (53) 



For the lapse with zero gradient at the puncture con- 
structed as a superposition of the odd and even lapse, 
the lapse is again one at infinity in order to measure 
proper time. The niultiplicator function $ appearing in 
the linear combination H27|) can be expanded in S as 



q-*0 



2x-M 
3 AP 



6^ + 0(6^) 



0.4714— +0(<53). 



(54) 



Therefore, at leading order the puncture lapse at late 
times is the average of the odd and even lapse. Next one 
can derive statements for the puncture lapse from "left 
to right" in the direction of increasing lapse. It turns out 
that the "zgp" lapse collapses at the puncture, 



the prefactor of the leading term of order 0{S) has a value 
which is half the one found for even boundary conditions. 
Finally, for the right-hand event horizon we obtain 



G 



lim 



X- 



q-*0 



1 



16 288V2 



7^6 In 



-330 - 72\/6 - 12\/3 
J2 



0.3248 -h 0.3967 



2r± 



9V6 + 22 
3x/2 + 4 

;3^ 



OiS'), (58) 



where again the finite value ^=^^ is approached asymp- 

totically. The reader should compare these results with 
the corresponding ones obtained for even boundary con- 
ditions. 



E. Analytical and numerical results for the 
Schwarzschild puncture lapse 



a. 



\x=0 



X-M 
2V2 52 



3 

0.9428 



M2 
f_ 

AP 



+ 0(5') 
+ 0{5% 



and at the left-hand event horizon, 

2 



zgp \r+ 



X- 



^AI 2' 



1 



288\/2 V 
-f330 



( 7\/6 



In 



9^/6 + 22 
3\/2 + 4 

?2 



60x/3)A_+O(<5^) 



1.1359— H-0(<53) 



in order Oip'). At the throat, however, since the odd 
lapse vanishes there and does not contribute to the "late 
time average" , with 



'■zgp kc 



9^0 



X-(M-A^) 
0.4714— 4- 0(,52) 



0((52) 



The late time expansions for the even and the "zgp" 
lapse have been tested both analytically and numerically. 

In Fig. 12 analytic results for the lapse are shown cor- 
responding to different values of the charge parameter q 
and both even and "zgp" boundary conditions. Since the 
(55) even lapse at the puncture is identically one it has not 
been plotted. We want to point out that the analytic so- 
lutions a^e„, and aj^p, if^ . have been evaluated 
at the appropriate values of r in the manner described in 
[T^ to produce the lines up to the box symbols, where 
the algorithm breaks down. Note that these lines go over 
smoothly when making use of the appropriate late time 
expansions of Subsccs . IIII CI and IIII Dl which verifies the 
leading order terms of our LTA. 

In addition, we have performed numerical simulations 
with "zgp" maximal slicing and vanishing shift. The 
particular results presented here were obtained with the 
Cartoon method described in [T9l | and implemented for 
the present purpose for s phe rical symmetry in the three- 
dimensional BAM code [lOj (our previous implementa- 
tions were for axisymmetry). The elliptic maximal slic- 
ing equation is solved with the BiCGSTAB algorithm, 
and for the details of the puncture evolution method see 
0. There are 512 equidistant points in the radial direc- 
tion X with a grid spacing of Aa; = 0.025i\/ and the outer 
boundary at about x = 12. 8A'/. 

Since the puncture is staggered between grid points, 
we obtain the value of the lapse at x = by polynomial 
interpolation using two grid points on each side. In order 
to calculate the lapse by interpolation at the throat and 
the left-hand and right-hand event horizon, we first have 

(57) 

to locate those on the grid as a function of time, c.f. 



(56) 
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Fig. 5 of • For this task wc compute the Schwarzschild 
radial coordinate r from the pre factor of the angular part 
of the metric as e.g. in 0,|23|, and identify the throat 
xc on a certain slice as the minimum of r and the left- 
hand and right-hand event horizon as isosurfaces with 
r = Or+ = 2M. 

As one can infer from Fig. [21 the numerical results 
agree quite well up to Tzgp ~ 40iV/ with the analytically 
derived "zgp" Schwarzschild lapse and its late time ex- 
pansions at the puncture and at the left-hand event hori- 
zon. Later on deviations occur since by then both round- 
off errors and noise introduced by the outer boundary 
contribute significantly. 

The agreement at the throat and at the right-hand 
event horizon is less convincing, note the different time 
scale of the plots. But as the maximal slices approach at 
late times the limiting slice r = "^fCHm = ^T' asymptoti- 
cally, it is very difficult to determine the location xc of 
the throat as a minimum of the Schwarzschild radial coor- 
dinate r since the minimum is very shallow, c.f. the lower 
plot of Fig. 13 Furthermore, since the lapse rises quickly 
from zero to unity near the right-hand event horizon, c.f. 
the upper plot of Fig.|31 small deviations from the correct 
Xq result in significant error for the value of the lapse. 
Those errors dominate at later times when the difference 
between the lapse and its finite asymptotic value is mea- 
sured, which is only of order 0((5^) and therefore poses a 
stringent test. 



IV. SLICE STRETCHING FOR VANISHING 
SHIFT 

The foliation with maximal slices - shown in Kruskal- 
Szekeres and Carter-Penrose dia gra ms in for the 
Schwarzschild spacetime and in |l2l | for the Reissner- 
Nordstrom spacetime - is uniquely determined by the 
lapse function. However, wc arc still free to choose the 
shift function and adopt different spatial coordinates on 
the maximal slices. In order to distinguish the 4-metric 
in those spatial coordinates from the radial gauge, J^l, we 
denote in capital letters with B the shift and with Gij and 
G the 3-metric and its radial component, respectively. 

Demanding a vanishing shift throughout the evolution, 
B = 0, is a trivial choice of the shift function which has 
been used in numerical relativity for its simplicity and 
which we will study in the following. In the context of 
maximal slicing, K = 0, this implies that the determinant 
of the 3- metric is time-independent. Hence the singular- 
ity avoiding property comes to light as the variation of 
the local volume remains fixed |2ll |. 

We will use the time-independency of det{G,„„} to 
study in Subsecs. llV Al and B the volume radius coordi- 
nate p. As in we will restrict ourselves to even bound- 
ary conditions. Here we are able to verify at the event 
horizon a conjecture regarding the late time behavior of 
the lapse at "large" p formulated for Schwarzschild in 
this reference. 



Again in the context of B = 0, in Subsecs. IIV CI to 
F we look at radial isotropic grid coordinates x. Those 
have been constructed in [Tl| such that the 4-metric coin- 
cides at all times with output from a numerical evolution 
of puncture initial data. Using "from left to right" the 
puncture, the left-hand event horizon, the throat and the 
right-hand event horizon as markers, we are able to point 
out slice stretching effects in the limit of late times. 



A. Volume radius coordinates 

Vanishing shift implies Eulerian observers obtained by 
choosing coordinates p{C{T),r) such that the level sur- 
faces p = const are timelike cylinders orthogonal to the 
slicing. As pointed out in Q for the Schwarzschild space- 
time and generalized here in a straightforward way for 
the Reissner-Nordstrom spacetime, one of these timelike 
cylinders, namely the one given by r = rc, is already 
known. By using the following expression for the volume 
radius. 



p'(C(Te„e„),r) =3 



Vpciy) ' 



(59) 



on each slice such a coordinate is found for even boundary 
conditions. To see this one may as a preparatory work 
verify that the partial derivatives of p with respect to r 
and C are given by 



d 1 d f y'^ dy 1 



9r dr J ^PcM ^/pc{r) 



(60) 



and 



a _ 1 9 [ yUy C d 



respectively. In deriving H51|l the reader should remember 
the notation teven{C, r) = Hc{r), l{TH|l. Furthermore, the 
expression 



oo 4 



y dy 

pc{y)^ 



the formula 



drc C 

dC --2rUrc-lM + ^) 



(62) 



(63) 



obtained by differentiating pc{i'c) = with respect to C 
and the statement (|21|l are helpful. 

Hence for even boundary conditions the Reissner- 
Nordstrom metric can be written as 



ds^ 



- ry rJ'T- 



even ^ even 



(^^y dp^+r^dn\ (64) 



9 



This result stated in Q one can immediately check re- 
membering the formulas Hl()|) . © and H2U|) for lapse, shift 
and time at infinity in the even case. One should observe 
that in 164|l the determinant of the 3-metric is given inde- 
pendently of r{C{Tf,ven),x) by det {Gmn} 

hence is time-independent. 



p"* sin^ 9 and 



B. The location of the event horizon in the limit of 
late times 

Next we want to show that a conjecture formulated by 
R. Beig and N. O Murchadha in ^ holds at the event 
horizon. For the Schwarzschild spacetime the conjecture 
was that the collapse of the even lapse should take place 
according to 



°ae«en(Te«en,°p) = -B("p) CXp 



on 



+C'(B2("p)cxp 
where p) behaves for large like 

M3 



Bi^p) ~ const • cosh 



qO/7,. on 



\ 3' 



M 



), (65) 



(66) 



Here the constant appearing in the equation for p) 

is given by 3°CHm°f^et)era/-M'^ = fgv^j hence note a mis- 
print in the previously stated reference |22l |. The partic- 
ular form of B{^ p) in (|66|l is motivated by the solution 
to the lapse equation, Aa = Ra, on the limiting slice 

We now want to look at the event horizon and first 
work out the late time behavior of the volume radius 
there to be denoted by ^ p^, the subscript C marking its 
time dependence. A good starting point for this task is 
Eq. which using the chain rule, ff = can be 

written as 



9V 

d5 



C 



'Pi' 



'Pi' 



0„2 



2(°r+ 



|M) 



ic^c(V) f (67) 



using (|21|) evaluated at °r+. To solve (in leading or- 
der) this ordinary differential equation (ODE), one can 
make use of (5-expansions for Kc{+), (|158|l with (|174|l . C, 
(I121|l . and hence (|122|) . derived in the Appendix. It 
is important to keep in mind that S at late times is decay- 
ing exponentially with time according to H46|) . With the 
first term in the second set of large parentheses of H67|) 
being of order 0(1) only whereas i^c(°r+) = 0{j^), it 
turns out that 



d"pt' 
dS 



1 



0{1) 



(68) 



is found with 
leading order 



\V3AP, (HH). We observe that in 



'pU^) ^ i/-3"Cu,n''^even In [S] as (5 (69) 
and hence 

' oo (70) 



' PciTeven) — i/S^Cl 



ira ' even 



as Tg 



hold. So indeed the conjecture of Q can be applied for 
the event horizon as the latter is found at large values of 
at late times. This is a slice stretching effect as will 
be pointed out for similar coordinates in Subsecs. IIV Dl 
and E. 

Now putting this result into the conjecture H65|l . we 
find in leading order 



^even \ ' even 7 



p+(Te„e„)) ~ 0(1), 



(71) 



SO the lapse approaches a finite value there in the limit 
Teven ~^ OO. This clcarly is in agreement with H5U|) . which 
states that the lapse at the right-hand event horizon has 
to approach the finite value jqV^ for maximal slicing of 
Schwarzschild with even boundary conditions. 



C. Isotropic grid coordinates 

In order to compare with numerical output calculated 
for maximal slicing of puncture evolutions, we construct 
isotropic grid coordinates as in ^J]. To focus on the 
dynamical features of the 3-metric rather than the static 
singularity at the puncture, we shall write 



'^^^ds^ =Gdx^ 



dn^ = 



g dx^ ^ J dfi^ 



, (72) 



introducing with small letters the rescaled 3-metric gij 
and its radial component g by 



and 



9{C{t),x)^ 



G{C{t),x) 
*4(a:) ' 



(73) 



(74) 



respectively. Here as pointed out in the conformal 
factor is determined on the initial slice from the 3-metric 
written in isotropic coordinates. Hence '^(x) is constant 
in time and given by H13|l . 



*4(a-) : 

So initially 

G{C = T = Q,x) 



mV qA'P 
2x1 4x2 



(75) 



^'^{x) and g{C = t = 0, x) = 1 



(76) 
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and in addition 



we obtain 



Q,x) = x^'^{x) 



(77) 



has to hold since on the initial slice the radial grid co- 
ordinate X is identical to the Schwarzschild isotropic co- 
ordinate R defined in (|ll|l in terms of the Schwarzschild 
radial coordinate r. 

We want to emphasize that when odd boundary con- 
ditions are used in a numerical simulation, the isotropic 
lapse (|12|) together with the formulas (|76|) and (|77|l are 
valid at all times and the code has to reproduce the static 
Reissner-Nordstrom metric in isotropic coordinates. 

For other boundary conditions, however, the compo- 
nents of the 4-metric evolve with time. Using that for 
vanishing shift the determinant of the 3-metric has to be 
time- independent, by imposing the initial conditions l|76() 
and lf77|l . in we were able to find 



9{C{r) 



x^-^^{x) 
KC{t),x 



(78) 



It is from G = ^-^g = 3^ = 7 then trivial to in 



fer for fixed slice label C the ODE 



dr 
dx 



(79) 



relating the Schwarzschild radial coordinate r to the ra- 
dial grid coordinate x. This equation can be integrated 
using the throat as lower integration limit by 



Vpc{y) 



= ± / 2/2*^(2/) dy. 



(80) 



where xc denotes the location of the throat on the grid 
as a function of time and the or "-" sign applies 
for the right-hand or left-hand side of the throat, respec- 
tively. One may readily check that for C = the grid 
coordinates coincide with the Schwarzschild isotropic co- 
ordinates where the throat is found at For later 
times, though, xc has to be found from H80|l by demand- 
ing that the coordinate transformation r = r{C{T), x) is 
consistent with the requirement of a vanishing shift. The 
latter can be written as 



B 



^ ^ dr\ dr _ ^ 
^ dr } dx 



(81) 



Remembering the formula (jSJ for (3 and writing the lapse 



as a 



y/pcir) gt dC 



dC dr 

dr 



, we observe that 



r4 dC dr 



(82) 



has to hold. By differentiating (|80|l with respect to C, 
which making use of H21() . H62|l and H63|l yields 



dr 



C 



2(r-|Af+^)V^ 



-Kc{r) 



Txh^'i^c)^, (83) 



dr 
d^ 



^2{r- 



1 



dxc \ dC 



-Kc{r)±-x^c'^\xc)-j^) — 



dC J dr 



(84) 



as a further equation for We want to point out that 
since in (|82|l the lapse enters, the behavior of the throat, 
to be determined by a comparison of 182|) and H84|l . de- 
pends on the boundary conditions imposed on the lapse. 
In Subsecs. IIV Dl and IIV El we show this concept for even 
and "zgp" boundary conditions explicitly when analyz- 
ing slice stretching effects in the limit of late times. A 
study of slice stretching in the context of more general 
boundary conditions is the subject of [23j . 



D. Slice stretching effects for even boundary 
conditions 



1. Throat xc 

Comparing (|82|l and H84|l . for even boundary condi- 
tions with the lapse given by IjlSfl or the partial deriva- 
tive of the height function with respect to C by (|21|l . we 
observe that 



dxc 
dC 







(85) 



holds and the throat hence remains at xc=o — x+. Re- 
membering that S for even boundary conditions has been 
found to decay exponentially with Teven on the funda- 
mental time scale il, (|46f) . we observe that the 3- metric 
"freezes" quickly at the throat and with H78|) for the 
rescaled radial metric component 



0< lim g {Teven, X+) 



and for the angular part 



_*8(a;+) 



< oo, 



< lim 



r2(T, 



even i ) 



< OO 



(86) 



(87) 



arc found. In particular, we want to point out that for 
the Schwarzschild spacctime g at the throat only grows 
up to the finite value (!)■*« 3.1605. 



2. Event horizon at xj, and x^ 



When discussing the location of the event horizon, note 
that with the throat remaining at a;-|_ for all times, in l|80|) 
the isometry (|14ll is present which maps the left-hand side 
of the throat to the right-hand side and vice versa. 

In order to derive an ODE for the location of the right- 
hand event horizon, we start by differentiating H80|) with 
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the integration limit given by the right-hand event hori- 
zon, i.e. the pair {r_|_,a;J}. By a calculation similar to 
the one leading to 1)61(1 or (jHSJ, it follows that 



i^CKc{r+) 



2(r+ 



2r+ 



~ c ^'^^> dC ■ 



Using the chain rule 



dC 



^ , and (5-expansions for 
Kcir+), with fT7^ . C, (Cni, and hence ^, ((n^ . 

derived in the Appendix, we find in leading order of 6 the 
ODE 



1 



0{l) 



(89) 



Note the similarity to the ODE (|68() obtained for volume 
radius coordinates. There we found that the event hori- 
zon moves to the right toward infinity, and in this limit 
the conformal factor approaches one. Presupposing this 
result and setting the conformal factor appearing in 189|) 
to unity. 



1 



, 2 dX(j 



d5 



(90) 



we infer that the right-hand event horizon expanded in 6 
as 



remains finite, on the left-hand side the corresponding 
term, 



hm , , _ , = 0, 



(96) 



approaches zero as the conformal factor grows without 
bound. 



E. Slice stretching effects for zgp boundary 
conditions 

Making use of the puncture data line element, H72() 
together with H78|) . and the coordinate transformation 
relating the radial gauge to isotropic grid coordinates, 
in [H we were able to verify the analytically derived 
4-metric with numerical output for maximal slicing of a 
Schwarzschild black hole with "zgp" boundary conditions 
and vanishing shift. In Fig. Owe show a^gp, g (plotted on 
a logarithmic scale) and r obtained numerically at times 
T^gp = {0, IM, lOM, lOOM}. The reader should observe 
that with the throat and the left-hand and right-hand 
event horizon acting as markers on the profiles of the 
3-metric (the "open boxes" and "triangles" in Fig. (JJ, 
by determining their locations on the grid as a function 
of time and evaluating the 3- metric at those locations, 
one can describe in the late time limit the effects of slice 
stretching. 



x;t((5) ~ ^-3CH,n^\n[5] as S~^0 
moves at late times to infinite values of x like 





•^n i'^even ) 



tm ' even 



cLS Tg' 



(91) 



(92) 



Note that the right-hand event horizon is found on the 
grid at a location that in leading order is proportional to 
the third root of the time. 

Furthermore, making use of the isometry, the left-hand 
event horizon is found to move toward the puncture ac- 
cording to 



(1 - q)A'P 



(93) 



zm ' even 



1. Throat xc 

For "zgp" boundary conditions, the ODE which deter- 
mines the location of the throat is most easily derived by 
comparing the partial derivative with respect to C of the 
"zgp" height function, H3U|) . appearing in H82|l with the 
term in the brackets of H84() . This yields 

-\c{Kc{^) + A) = xl^\xc)^. (97) 

Therefore, for xc{5) based on iS-cxpansions for Kc{oo), 
(CSHll, C, ((nT|l . and hence ^, ifT^ . in leading order of 
5 the ODE 



For g at both the left-hand and right-hand event horizon 
in leading order the statement 



-Curn^ \ + 0{l) = xl-^\xc)^ (98) 



9{' 



even-} -^(j ('^euen )) 



{j^Clira'^even^ ^ 



0jS Tp. 



(94) 

can be made. Hence at the event horizon the rescalcd 
radial metric component is proportional to Ttien and 
grows without bounds. For the angular part we note 
that whereas on the right-hand side the expression 



is obtained. Note the similarity of this ODE for the 
outward-movement of the throat in isotropic grid coor- 
dinates at late times to the ODE 168|l describing the lo- 
cation of the event horizon in volume radius coordinates 
and to H89|) for the right-hand event horizon in the con- 
text of even boundary conditions. Neglecting again the 
conformal factor in order to infer the behavior of xc(5) 
in the limit of late times, by solving 
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wc find 



that the right-hand event horizon expanded in 5 as 



xciS) ~ ^-3C;.„„171n[(5] as 5^0. 



(100) 



+ {S) ~ ^-6Ci„n^ln[S] as (5 ^ 



(107) 



Since 5 is decaying exponentiaUy with r^gp according to 
(I51|l on the time scale ^Izgp = 251, H52|l . the location of 
the throat on the grid at late times is hence dominated 

by 



2 ^lim^zgp '^zgp 



(101) 



Of course, this simple solution is not expected to give a 
good approximation at early times since both the expan- 
sion in 5 and neglecting the term ^'^(ccc) are not valid in 
this regime. But at late times the numerically computed 
outward movement of the throat agrees well with the 
power law (|101|l obtained for Schwarzschild and "zgp" 
boundary conditions, c.f. Fig. ^ 

Since we arc interested in the late time behavior of the 
3- metric, bearing in mind the definition of 5, (|36() . we can 
infer from (|78|l that the rescaled radial component of the 

metric at the throat blows up as Tzgp, i.e. 



g{Tzgp,Xc{Tzgp)) 



\ 2 ^^'^'^ 



as Tzgp — > oo. (102) 



The angular part, though, remains finite. 



< hm 



r'^{Tzgp,Xc{Tzgp)) 



r.gp^oo ^^{Xcijzgp)) 



rl^^^^ < oo. (103) 



2. Right-hand event horizon 



Now differentiating (|80|l with the integration limit be- 
ing the right-hand event horizon, i.e. the pair {r+,x'^}, 
with respect to C, it follows that 



moves at late times to infinite values of x like 



X'fji'Tzgp) — \/ 3Cli 



as Tzgp > oo. 



(108) 



i.e. with a location on the grid again proportional to the 
third root of the time. However looking at the "collapse 
speed" , i.e. the rate at which the numerical grid falls into 
the black hole, one can see from 



lim ^ ''7'"' ^ 



(109) 



that the right-hand event horizon moves at late times by 
a factor of « 1.2599 faster than the throat. For the 
3-metric, 



9{'^zgp 1 ^'ci'^zgv)) 
and 



{^^ClimTzgp^ ' 



zgp OO (110) 



f {Tzgp, Xq {Tzgp)) 



0< lim • v^gp.-^cV-.gp; 

r.„-.oo ■^^x+{Tzgp)) 



r I < OO 



(111) 



can be found. So the imboimded growth of g is again 

proportional 
Calculating 



proportional to Tzgp and the radial part remains finite 



dg{Tzgp,X^,{T^gp)) 
rl-r 

lim 



(112) 



we claim that for Schwarzschild the blowup of g should 

occur for the right-hand event horizon by a factor 
3^4 



23 



0.7973 slower than for the throat. 



1 



2{r, 



^CKc{r+) 



. dXr 



dxc 



(104) 



Adding the equations H97|) and (|104|l and using 
Kc{oo) ~ Kc{r+) =ri + 0{6^), itTT^ . we find in leading 
order of S the ODE 



-2Q 



lim 



(105) 



Note that this ODE for the location of the right-hand 
event horizon at late times only differs by a factor of two 
from H98|l obtained for the throat. As for the throat we 
set the conformal factor to one and infer from 



1 



, 2 dXQ 



dS 



(106) 



3. Left-hand event horizon 

For the left-hand event horizon one has to take care of 
the minus sign on the right-hand side of H80|) . As for the 
right-hand event horizon, with 

rl 1 



2(r+-|M+^) 



7;CKc{r+) 



- -^c'*'(-c)^+4*^(-c)^ (113) 
subtracted from H97|l we obtain 

( 1^ . 2 \ ^ 



+0{5^) = (114) 
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As the prefactor of 5 on the left-hand side of H114() is 
positive, we observe that the value of is decreasing 
as 5 is decaying to zero. For this reason wc do not set 
the conformal factor to unity here but observe that in 
the limit of late times the left-hand event horizon will be 
found at a finite location on the grid to the left of its 
initial location, i.e. 

0< lim a;p(r^gp) = <x+. (115) 



It then follows that the 3-metric "freezes" as 
< _ lim_ 5(T2gp,Xp(r^gp)) = — ^i^^^— j 



and 



< oo, 
(116) 



< lim ''^^^^av^^'ci^^av)) ^ ^ jl ^ ^ ^^^^^ 



holds. 



4-. Peak in the radial metric component at 

Finally we want to discuss the peak occurring in the 
radial component of the rescaled 3-metric at location x^. 
Inferring from Fig.|3|that Xq is found to lie in between the 
location of the throat and the right-hand event horizon, 



(118) 



we claim that upper and lower bounds for the location of 
the peak are given by (|101() and (|108|l . Furthermore, the 
corresponding value of the Schwarzschild radius is then 
bounded from below by rc,.^ and from above by r-|_. 

From H78|) we can therefore infer the following bounds 
for the peak of g, 



^4 



< 



(119) 

which - setting for simplicity the conformal factor to 1 - 
have to hold in leading order at late times. 



event horizon on a given slice had to be determined. 
Whereas an event horizon can be located by searching for 
the isosurfaces r = '^r+ = 2M, searching for the throat as 
the minimum of r is difficult at late times as the mini- 
mum is very shallow, c.f. Fig. |3| Instead the throat has 
been located via its value of the lapse function that is 
given at late times by (|57|l . Finding the event horizon by 
the corresponding values of the lapse yields locations that 
are very similar to the ones obtained with the isosurface 
based method. 

Note in Fig.^that for the left-hand event horizon both 
the location and the value of g are essentially constant in 
agreement with analytical predictions. In addition, one 
can infer for the throat and the right-hand event hori- 
zon that the late time behavior for the locations is very 
well reproduced by H1U1|I and (|108|1 . The correspond- 
ing curves have been fitted to the numerical value at 
Tzgp = 80M since the LTA only reproduces the behavior 
at late times but is not capable of describing the behavior 
at early times. Fitting at this value of r^gp is a compro- 
mise between applying the LTA at reasonably late times 
and avoiding large numerical errors at later times. 

For the blowup of g the accelerating growth taking 
place at the throat and the right-hand even horizon is re- 
produced qualitatively by (|102|) and H11U|) . However, the 
reader should bear in mind several factors which con- 
tribute to deviations. On the one hand several simpli- 
fications, valid only at late times, have been used when 
deriving the late time statements. On the other hand the 
numerical errors in determining values at only approx- 
imately known positions by interpolation in a sharply 
peaked metric profile are not negligible. Furthermore, 
the numerical outer boundary is an approximation to a 
radiative boundary condition, and its infiuence is felt by 
the interior solution at later times, producing an error in 
the height of the metric peak. 

Finally, Fig.^also shows the bounds for the peak in the 
radial metric component derived in Sec. IIV E 41 These 
bounds are not very sharp, c.f. the gray area in Fig. 21 
But note that the peak in the radial metric component 

has to grow proportional to Tzgp . An exponential growth 
of the 3-metric, which in numerical relativity was some- 
times expected for the slice stretching problem, can be 
ruled out. 



Numerical results 



CONCLUSION AND OUTLOOK 



The late time statements about slice stretching for van- 
ishing shift have been compared for "zgp" boundary con- 
ditions and vanishing charge to numerical simulations, 
see Sec, nil EI for comments on the numerical details. We 
summarize the result in Fig.0] Note that this data helps 
explain e.g. Fig. 3.13 for the outward movement of the 
event horizon and Fig. 3.16 for the growing peak in the 
radial metric function g in |24l |. 

To determine the numerical data points of Fig. ^ the 
locations of the throat and the left-hand and right-hand 



In this paper we have considered the late time limit 
of maximal slicing of the extended Reissner-Nordstrom 
spacetime for both even and numerically motivated "zgp" 
boundary conditions. In this limit, the leading order 
terms of the lapse function and the 3-metric at the punc- 
ture, the left-hand event horizon, the throat, and the 
right-hand event horizon together with the correspond- 
ing locations on the grid as a function of time have been 
obtained for vanishing shift. 

Taken together, the late time analysis for the lapse 
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and statements for the grid locations provide an analytic 
understanding of the numerically observed lapse profile 
for "zgp" boundary conditions with the "collapse of the 
lapse" and the "outward moving shoulder". The lapse 
has been shown to collapse at late times at the punc- 
ture, the left-hand event horizon and the throat, with 
the latter moving outward. The value of the lapse at 
the right-hand event horizon, however, moving outward 
slightly faster than the throat, approaches the finite value 
< %L < °.Ci^ ^ 3 ^ _ 0.3248. 

Furthermore, we have discussed slice stretching effects, 
which have unpleasant consequences in numerical simu- 
lations. The associated effects, such as outward-drifting 
of the throat and the right-hand event horizon (slice 
sucking) and unbounded growth of components of the 
3-metric (slice wrapping), have been described analyti- 
cally. In particular, we were able to rule out an exponen- 
tial growth of the latter. 

For the "zgp" Schwarzschild case, the analytic state- 
ments for the outward moving shoulder in the lapse pro- 
file and the development of a rapidly growing peak in the 
rescaled radial metric function have been verified numer- 
ically. 

Note that for constant mean curvature slicing of the 
Schwarzschild spacetimc the time scale of the collapse 
of the lapse has been determined as a function of K in 
|25j . It would be interesting to extend our studies in this 
direction. 

For more accurate and longer lasting simulations with 
singularity avoiding slicing, promising strategies are ei- 
ther to use singularity excision or to use a non-trivial 
shift, which as the slice stretching develops reacts by 
pulling out grid points from the inner region. This has to 
be done based on appropriate conditions imposed on the 
3-metric. Candidates like the minimal distortion shift 
have been studied numerically as early as e.g . IS l. The 
Gamma-freezing shift conditions discussed in |9l. \2(L 
have allowed evolution times of lOOOM and more, ap- 
proximately freezing slice stretching. It is therefore of 
significant interest to provide an analytic understanding 
of the late time behavior of these shift conditions as is 
now available for maximal slicing. Some analytical work 
has been already done in [T^ . and we plan to examine 
shift conditions in more detail in the future. 
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APPENDIX: 
PROOF OF THE LATE TIME ANALYSIS 

1. Ansatz 

For our LTA, in contrast to Q we will not rescale 
quantities like C, r or 5 by their units of M here, since 
quite often the involved powers of M give insight into 
the underlying physics and provide a first check of the 
formulas. So throughout this proof M will be under- 
stood as a positive constant not necessarily equal to one, 
M > 0. Furthermore, we will restrict the charge param- 
etcr q ~ the squared ratio between charge Q and 
mass M , to < q < I to discuss only the non-extremal 
Reissner-Nordstrom spacetime. 

The ansatz for the LTA is based on an expansion in 

5 = rc — fcii^, (ESI) I i-e. the difference between the ra- 
dial coordinate at the throat rc on a slice labeled with 
C and its limiting value rc,.^^ corresponding to Cum- As 

6 for the non-extremal Reissner-Nordstrom metric de- 
creases strictly monotonically for C € [0, Cum], the range 
of 5 is (5 S [0,tT+], and the late time limit is obtained by 
considering the limit S ^ 0. Here for convenience Ci <0 
as in (|37|) has been introduced. Note that cr+ as a func- 
tion of q decreases from ^ for the Schwarzschild case 
to zero for the extremal Reissner-Nordstrom case, where 
the LTA obviously breaks down. 

We want to point out that, for the non-extremal 
Reissner-Nordstrom metric, S defined by (|36|l yields a 
valid ansatz for the discussion of the late time behavior 
of spacelike maximal slices extending on both sides of the 
throat to infinity (referred to as "horizon-horizon" slices 
in unless one starts with the initial slice being the 
limiting slice r = rc^^ ■ 

The task of finding the (5-series for the times at infin- 
ity and the lapse at the puncture (a~ |2;=o)j the throat 
(a* |rc)i ^-iid the event horizon {a^ \r^) can now be re- 
duced to expanding C and the integrals Hc{oo), Kc{r+), 
and Kc{oo) in 5. This one can see from the formulas for 
the r's and a's, ^ and for the even 

and "zgp" boundary conditions. 

2. Expanding the slice label C{5) 

With rc the unique (double counting for C = Cum) 
root of the polynomial pc{r), it follows that 



C = rcy-rl, + 2Mrc - qAP 
= rcy^r+ - rc^Jrc - r_ 
= {^ + rc,^J^Jo^-8^b-G-. (120) 

Here has been used and without loss of generality 
the positive root has been chosen to yield < C < Cum , 
since the negative root, —Cum < C < 0, together with an 
inversion of the time direction, r — > — r, would result in 
the same maximal foliation. Making use of the positive 
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constant ^ defined in (|39|l and expanding C up to second 
order in S yields 



C = Cii 



2^-cr+cr_ 



S' + OiS-"). (121) 



Here Cum — rc^ .^-^— cr+cr_ can be read off from H120() in 
the limit S ^ 0. For later use (|121|l is diff'erentiated with 
respect to 5 to obtain 



dC M - /i 



(5 + 0((52), 



an upper bound for Cum — C is given by 



< Cum ~ C < 



lim c-2 _ Clirii ^2 



\c=o 



a. 



2 " ' 



(122) 



(123) 



and also the difference between C^^j^ and is calculated, 

Cl^^C' ^S\5'+^i5 + ,.), (124) 
using as in (^UJ- 

3. Expanding the time at infinity t{S) 
a. First Task 

In order to prove the statements for the exponential 
decay of d with time at infinity, one has to show for even 
boundary conditions that 



Hc{oo) = -nin[S]+ A + 0{d) 



(125) 



holds, which inverted with the correct expressions for 
and A as in H38|) and H45|) implies (|46|) . Once this is shown 
for the even case, together with (|120|) the result for the 
"zgp" case follows trivially from (|32|l since 



r.gp = 2Hc{^)-— = -2n\ii[S] + i2A- ^) + 0{S). 

(126) 

The basic feature of the proof now consists of finding ap- 
proximations Ti for Hc{oo) which, being straightforward 
to integrate and to expand in S, allow estimates for the 
remainder in order to complete the task of showing (|125|l . 
The necessary preparations and estimates will be derived 
in the following. 



b. Preparations 



We start with 

Fc(oo) 



C{6) ds 



f{r)\/pc{r) 



(127) 



obtained by substituting 

r = s + rcii„ 



(128) 



in the integral p9|l to be calculated for the time at infin- 
ity if^ . and write C{S) as in H120() . For the term 

one can then find with A± <0 as in (|41|l the following 
expansion into a partial fraction. 



fir) 



(s + rc„J'^ - 2M (s + rc:,J + qAP 
A, A_ 



1 



(129) 



Also one should observe that for the polynomial pcii,„ {r) 
corresponding to the limiting value Cum with rc, .^^ being 
a double root. 
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lim 



= {r~rc,,y{r^ + 2irc„ 
one can by substituting (|128() write 



(130) 



M)r - 



(^2 



(131) 



Hence together with (|124|) we obtain the expression 

Pc{r) = Pc,.Jr)-Cf,„,+C' 

= s^{s^ + fis + ly) ~ 6^{d^ + fiS + ly). (132) 



With 1.1, u > the inequalities 
< iy{s^ - 6^) < pc{r) 



(133) 



pcir) < [s^ ^S''){iy + 2{s^ + fis)) (134) 

hold for < S < s, which allow us to derive the two basic 
estimates needed in what follows. The first. 



1 



VpRO V^(s2 - 52) 



< 



< 



< 



y/iy{s'^ - d'^)pc{r) 
y/ V + 2{s'^ + lis) - ^ 
^/vpc{r) 



l + f(s2 + ;xs)-l 
^i^(s2 - 5^) 



(52 



(135) 



is found by applying (| 134(1 in the numerator, ((133() in the 
denominator and using < 1 -|- f , a; > 0. With 



d_ 
dC 



C 



\Jpc{r) 



pcir)i 



(136) 
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the mean value theorem can be apphed to yield together 
with H123|) the inequality 



late time limit 5 — > as 



C 



fir) \VPc{r) VpcZJr) 



< 



Cii 



/x2 / \ a ' 



(137) 

which by continuity is also valid for ?■ = r+. Using (|128|l . 
and also this can be written in the form 



1 / C'iS) 



Cm 



fir) V ^/pc{r) VpcZJr) 



< 



< 



PG(rY- 



(138) to yield 



ds 



(143) 



--ln[<5]+ln [2^]+0{6), 



I {5a" 4 



(s — a±)\/ ~ 



ds 



(144) 



It is worth mentioning that for the two estimates p35|l 
and contact with Eq. (3.12) and (3.13) in 13 can 

be made by recalling the values "/i = 4M and = 
for zero charge. 



Cli 



l-(^ + — 

cr+ (J- 



ln[6] 
In 



(145) 



0(1). 



c. Estimate for the integral Hc{<x) 

The integration domain for the integration taking place 
in 1127(1 is split into A : < S < s < y^Sa+ < cr+ and 
B ; \/oaZ < s < oo. So we can write 



Similarly, for TC^ written as 



when introducing ^± > as in lO, with 



(146) 



introducing the approximation integrals Ti^ and Ti^ 
given by 



Cum ds 



f{T)^y{s^-5^y 



Clim ds 



/ (5(J4 



(140) 
(141) 



ds 



rB± 



.J_ln[5]-^ln 



ds 



2y/v + /i 



{s-cf±)^Jpcu^{r) 



ds 



1 



2^o-± 



A. 



ln[<5] 
o(l), 



In 



207 + ^ 



(147) 

o(l), 
(148) 



4z^ 



Using H129() we obtain H-^ as 



we obtain 



(142) 



where the following inte gral s can be computed explicitly 
(using e.g. Chap. 21 of |23) and expanded in S in the 



(149) 



l-(^ + -))ln 



-Cum{i++^-)+o{l). 



2v^ + /i 



41/ 
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Estimating the remainders for AH^ 
yields 



iJ^ -7^^ now 



AH' 



1 



C{8) 



Cu 



\ds 



< 



1 C{8)-Cu, 
/W ^Jvc{r) 



ds 



1 5(7^ 



Cu 



< const 



: ds 



= 0{V6). 



Ids 



(150) 



Here (|121|) has been used to neglect the first and H135() 
to bound the second summand. Given the second basic 
estimate AH^ = \H^ -H^] in the limit 6^0 

has a bound of the form 



1 



C{6) 



Clir. 



d. Second task 

Next, we have to show that in the limit the last 

equality in 



dHc(oo) 1 , ,dC 1 , , , , 



holds. This expression is obtained by making contact 
between He {oo) and Kc (oo) using (|^ . Note that inte- 
grating (|156|l . which yields 



Hcioo) = -nin[6]+A*+0{S), 



(157) 



and inferring A = A* finally proves the original task for- 
mulated in 1123). So with itT^ one has to show 



(158) 



e. Estimate for the integrals Kc{r+) and Kc{oo) 



< const • S 



(s -I- rci„J^ ds 
(5^(5^ + lis + 



(151) 



Substituting (|128|1 in ()17|l . we can write Kc{oo) as 



Following and splitting this integral into 



(5 + rci„„Y ds 

(s2(s2 ^ ^5 ^ jy)); 



< const • 



(.S2 - ^2)1 



(152) 



and 



(s -I- rc,.^)^ ds 



< 



00 

/■ ds 

const • / ^ = C'(l), (153) 



we finally find 



H^ -n' 



0{6). 



(154) 



Putting the results for the Ti's and the estimates of the 
remainders AH = \H - 1-L\, ((T151 [T^ and (jTCUl [T31jl . 

together we obtain 

^^0(00) = -01n[(5] + A + o(l) (155) 

with and A as proposed in (|38|l and (|45|1 . However, 
comparing (|155|l with the original task as formulated in 
it turns out that our estimates so far are (as in 
I ) "not quite good enough" to conclude the form of the 
higher order term appearing in the exponential decay of 5 
with r. The estimates give and A, but we now proceed 
to improve the estimate for the o(l) term. 



Xc(oo) 



X+ , X- \ ds 



making use of a partial fraction with x± defined in \A'6\ . 
For Kc{oo) the integration domain can be split at the 
event horizon, which is appropriate since one is later on 
also interested in Kc{r+). Note that this was not pos- 
sible for Heir) since this integral diverges at r+. Hence 
A : < S < s < a+ and B : (7+ < s < 00 are the integra- 
tion domains. Again we can write 

Kc (oo) = K^+K^ 

= [K"^ - K.^) + K.^ + (K^ -K.^) + K.^ (160) 

and introduce approximation integrals K,^ and /C^ by 



x+ 



ds 



(S - ^)2 S- 



00 



X- 



(161) 

ds 



l^s-^V s- 



J^'^ + x+J^'^ + x-J' 



(162) 
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Explicitly computing and expanding the J^s yields 

ds 



A* 



J 



jA+ ^ 



ds 



O(lnH), 



6 



(S-^)2^1/(S2_J2) 



= 0(ln [(?]), (164) 



J 



A- 



2 

ds 



= ^1+0(1). (165) 



Abbreviating In [• • • ] = In 
also 



cr+/J+2(l/+yi7^<T5j_+(T+/j+!/) 
a+(/^+2v^) 



ln[...], 



(s-^)^v^^:~F) 



(166) 
(167) 



16(2cr_ - M + 2y'(T2 ^ ^ 
(2cr_ - A'/)2(A/-//)2 
^ ^ ln[...], 



(M - ^)2 
ds 



(168) 



1 i^^ - Av 



are found. Summarizing results, 
^A X- 1 



ln[. 



(169) 



and with 77 defined already in (|44() we obtain 

Last, but not least, the remainders in the limit S ^ can 
be estimated for AK"^ = IK"^ - IC^\, 



(S - ^)2 S2 



1 



1 



ds 



Vp^^ v/^(.s2-52)^ 

< const. 7^|£==0(i), (171) 



using 1123 and for AK^ = |A'^ - /C-^I, 



(163) AK^ 



X+ 



is 



M^\2 
2 ) 



ds 



< const • 5 



^ f{r){s + rc,,jUs 

(s2(s2 + fiS + 



0((52),(172) 



using p38|l while considering (|121|l . Now taking finally 
into account the results for the /C's and the estimates 
AK =\K -IC\, (UnnrnZni) and ((miini, together with 



M - [I 



(173) 



the statements (|158|) and therefore H125() are proven. Fur- 
thermore, since A'c(oo) and Kc(t+) differ by only, 
we can conclude from (|170|l and (|172|l that 



(174) 



holds. We can use this expression to extend the study 
of since (|174|l allows us to obtain expansions in (5 not 
only at the throat but also at the event horizon. 



4. Expanding the multiplicator function $(5) and 
the lapse aiS) 

Using H158|l the multiplicator function $ can be ex- 
panded in 5 as 



2 Kcico) + 2 2x-M 



So at late times the "zgp" lapse is in leading order the 
average of the odd and even lapse. 

The lapse at the puncture is minus and plus one for 
the odd and even boundary conditions, respectively, and 
therefore at a; = the "zgp" lapse constructed as a linear 
combination is obtained as 

O^'zgv U=0 = * • a^en U=0 +(1 - ' OL~dd U=0 



= 2$ - 1 = -- 



X-M 



b'^^Oif). (176) 



Next, the value of the lapse at the throat tq is found 
for the even case from (|16|) using 136|l . 



a; 



1 (5 + rc,. 



even KG 



6^0(f), (177) 
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and for the "zgp" boundary conditions it is at leading 
order half of this value, 



'^zgp \rc ^ ' '^even \rc 



S + 0{S^), (178) 



since the odd lapse vanishes at the throat and docs not 
contribute in the late time average. 

Finally, the lapse at the event horizon can be expanded 
in S. Remembering that r+ is a root of /(r), for the odd 



lapse it follows that a^^^ |r+ = 
r_)_ is given by 



±-3-. The even lapse at 




) 6' + OiS^), (179) 



where the expansions (|12l)|) . H158|) . and (|174|l have been 
used. With the series (|175|) for <i> it turns out that for 
the "zgp" boundary conditions on the left-hand side of 
the throat the value of the lapse at the event horizon 
collapses to zero as 



*zgp 1'"+ 



$ • a 




S^ + OiS^). (180) 



On the right-hand side, though, with 




odd !'■+ 



)]6' + Oi6') (181) 

as for the even boundary conditions, the value "^'i"' is 

approached. The latter is an important example for the 
conjecture regarding the late time limit of the lapse at 
the right-hand event horizon stated in our previous paper 

5. Comments 

We want to make several comments regarding the proof 
of the late time analysis. 

Since the whole proof depends on the statements 
(tn5|l and fT7^ for i?c(oo) and A'c(oo) - ifc(r+), 



these two fundamental expressions have been checked 
numerically in |l2l |. This can be done by calcu- 
lating for different values of q the integrals Hc{oo), 
Kc{r+), and Kc {00) with Mathematica and by 
verifying that both ^ \Hcioo) - {-D,ln[S] + A)\ and 
j2 \Kc {00) — Kc(r+) — 77 1 behave like 0(1) in the limit 
S ^ 0. Hence not only the formulas (I2H1 ESJ and 
found for fi, A and 77 as functions of q can be confirmed, 
but also the order of the "rest term" appearing in H125|) 
and 1)1 74|) can be verified. 

Note that astrophysically relevant black holes are ex- 
pected to have < q I as it would be very difficult 
for any astrophysical body to achieve and/or maintain 
a charge to mass ratio = ^-jf of greater than 10^^® 
|30l |. since otherwise particles of opposite charge would 
be attracted selectively. By analyzing the terms depend- 
ing on q appearing in the S expansions, one can see that 
for moderate charges these expressions deviate only very 
little from their Schwarzschild values denoted by an up- 
per index "0" referring to zero charge. For quantitative 
details and plots of Q, the k's, and rj as a function of q 
see [13. 

However, in the limit q ^ I the expressions for Q and 
A diverge like 



.1 1 



-.M 



and 



2 

Furthermore, with C;,;™ vanishing like 



(182) 



(183) 



(184) 



the pre-exponcntial factors approach zero for q ^ 1 as 



q->l 



q^l 



(185) 



and the LTA obviously breaks down. This, however, is no 
surprise as for the extremal Reissner-Nordstrom space- 
time containing a naked singularity our ansatz for 6, H36|l . 
is not valid. 

Using previously stated expressions, it should be 
straightforward for the reader interested in the full 4- 
metric at the event horizon to expand also the radial 
metric component, (|SJ|, and for even or "zgp" boundary 
conditions the shift, ((Sj). But note that both 7 and f3 
diverge at the throat. This, however, is only a coordi- 
nate effect of the radial gauge. For numerical purposes 
it is therefore preferable to perform as in Subsecs. IIV Dl 
and IIV El a late time expansion of the 3-metric (or the 
shift) in suitable spatial coordinates such as the isotropic 
(or, in the case of a non-vanishing shift, isothermal) grid 
coordinates. 
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FIG. 2: For puncture, left-hand event horizon, throat and 
right-hand event horizon from top to bottom the value of the 
lapse as a function of time at infinity is shown. The lines 
- plotted in gray or black corresponding to even or "zgp" 
boundary conditions - represent analytic results obtained for 
q = (solid), 0.5 (dashed) and 0.99 (dot-dashed). In ad- 
dition numerical results obtained for vanishing charge and 
"zgp" boundary conditions are plotted as symbols. Those are 
found to agree very well with our analytic statements at the 
puncture (circles) and the left-hand event horizon (downward 
pointing triangles). The deviations occurring at the throat 
(boxes) and the right-hand event horizon (upward pointing 
triangles) are explained in the text. 
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FIG. 3: For "zgp" boundary conditions numerical results at 
times r,gp = {0, IM, lOM, lOOM} are shown to demonstrate 
the "collapse of the lapse" and "slice wrapping and suck- 
ing" taking place as the limiting slice r — is approached. 
The "symptoms" of slice stretching are the development of a 
rapidly growing peak in the radial metric function (note the 
logarithmic scale!) and the outward-drifting of coordinate lo- 
cations. To study the latter, one should observe the location 
of the throat, the peak in the radial component of the rescaled 
3-metric, and the right-hand event horizon denoted by boxes, 
stars, and upward pointing triangles, respectively. 




FIG. 4: Top: Slice sucking is studied numerically by plot- 
ting as a function of time the location on the grid of the 
throat (shown as boxes) and the right-hand event horizon 
(denoted by upward pointing triangles). In between the lat- 
ter two the peak in the radial metric component (plotted as 
stars) is found. The analytically derived late time behavior 
(fitted at Tzgp = 80M and predicting a growth proportional to 
T"i/p with known prefactors) agrees very well with numerical 
results. Bottom: Similarly slice wrapping can be investigated 
by evaluating the radial part of the 3-metric at those loca- 
tions. At late times again agreement with the analytically 

4/3 

derived blow-up proportional to r^gp is found. The gray area 
indicates the bounds obtained for the peak in g. 



